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RESONANCE, PARTICLE TRAPPING, AND LANDAU DAMPING IN 
FINITE AMPLITUDE OBLIQUELY PROPAGATING WAVES 

P. J. Palmadesso 
Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

The equations of motion for a particle in resonance with a small finite ampli- 
tude wave (oi, k) are solved approximately, using secularity free perturbation 
theory. The wave propagates at an arbitrary angle to a uniform background 
magnetic field, B , in an infinite collisionless plasma. The wave fields include 
a longitudinal electrostatic component and elliptically polarized transverse ; 
electric and magnetic components. The trajectories of trapped and resonant 
untrapped particles are described, for each of the possible wave-particle 
resonances defined by the condition k z V z - N q B Q /m, where N is an 
integer. These trajectories are used to construct an estimate of the nonlinear 
time dependent Landau damping rate of the wave. 
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I. INTRODUCTION 


The trapping of particles by finite amplitude narrowband waves leads to a 
number of important nonlinear effects. Particle trapping is responsible for the 
saturation of the Landau damping/amplification mechanism in the nonlinear 
regime . 2 Electrostatic trapped particle instabilities , 3 " 6 and the triggering 
of V.L.F. emissions by finite amplitude whistlers in the magnetosphere 7 ’ 8 
provide other examples of effects in which trapped particles play a major role. 

Explanations of these phenomena depend heavily on a knowledge of the tra- 
jectories of individual resonant particles, i.e., trapped and nearly trapped 
particles, in the fields of the primary wave. In all of the examples mentioned 
above, the primary wave is either an electrostatic wave or a circularly polarized 
transverse wave propagating parallel to a background magnetic field. The 
trajectories of particles in fields of these types are well known . 1,9 

In this paper we obtain relatively simple approximate analytic expressions 
which describe the motion of single particles in resonance with a wave of some- 
what more general type. We consider a small finite amplitude wave propagating 
at an arbitrary angle to a uniform background magnetic field in an infinite 
collisionless plasma. The wave has both longitudinal electrostatic and trans- 
verse electromagnetic components. 

An immediate application of the results of the trajectory analysis is possible. 
The approximate orbit expressions for resonant motion in an oblique wave have a 
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form which is similar in some important respects to the corresponding expres- 
sions for motion in an electrostatic wave or a parallel propagating transverse 
wave. By recognizing this, one is able to obtain in a simple way an estimate 
of the time dependent collisionless damping coefficient for a finite amplitude 
oblique wave. One finds that after some initial exponential damping the wave 
executes damped amplitude fluctuations. In the absence of instabilities, the wave 
amplitude would ultimately become constant. 

We begin in Sec. II by solving the equations of motion for a single non- 
relativistic particle, using secularity-free perturbation theory. 10 The constants 
of motion and bounce frequencies are found, and the time dependence of the 
particle velocity is determined. Finally, in Sec. Ill, an estimate of the time 
dependent damping coefficient is obtained. 

II. PARTICLE ORBITS 

Consider an electromagnetic wave (w, k) propagating at an angle a to a 
magnetic field (B Q ) in a plasma. If the wave amplitude is not too large, it is 
reasonable to assume that the propagation characteristics of the wave are pri- 
marily determined by the nonresonant particles which form the bulk of the 
plasma, and that these particles are describable in terms of the linearized 
theory. Given the additional assumption that y / co < < 1, where y is the damping 
rate of the wave, the wave fields may be written in the form 1 
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Bw <R’ t) = ^x B lsin^ + (e x x e k ) B 2 cos 0, 


( 1 ) 


B w (R- O = &K Ek sin 6 + C" /k ) Bw* £k- (2) 

where 

0 = fe- B - " t, 

& k = fc/ k - 

and the coordinate system is oriented so that 

£x = (fex So) /k B 0 sin a - 
£y = (Bo * ^x)/ B 0’ 

® z - Bo /B 0- 

Two types of wave-particle resonance are possible with obliquely propagating 
waves. Cyclotron resonance occurs when the rotating wave fields seen by a 
moving particle turn with a frequency close to a multiple of the particles gyra- 
tion frequency, fi 0 ; i.e., when 

O) - k V + N fl„ = 0 

z z 0 

where V z is the component of the lab frame particle velocity parallel to B 0 , and 
N is an integer. In the simplest case (N = ±1), this means that the average 
value of the angle between the wave electric field and the perpendicular 12 velocity 
vector of the gyrating particle changes very slowly. The power delivered to the 
particle, E W >V, then has a slowly varying component, and large amounts of energy 
are exchanged. Particles for which K E w > \ is bounded are said to be trapped. 
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Note that it is necessary to speak of the average value of $E , when describ- 
ing resonance. The plane in which particles gyrate, in zero order, differs from 
both the plane of polarization of g w , and the planes of constant phase. 

Electrostatic resonance is possible if <^>/k z < c, where w/k z is the speed 
along B 0 of a point of constant phase. In a frame moving with velocity (<Vk z ) e , 
the wave has only a longitudinal electrostatic component and a transverse mag- 
netic component, and both are stationary. Particles nearly at rest in this frame 
( a) - k z V z ~ 0) experience electrostatic forces which vary only slowly, so that 
E w - V again has a slowly varying component. Such particles may be trapped in 
the potential wells of the longitudinal field. 

Equations of Motion 
If we introduce the velocity variable 

v = V - (to/ k ) e (3) 

'■'o 'V. v z ' Z 

and use (1) and (2), the equation of motion 

m i = q + Vx (B 0 +g w )} 

can be written in the form 

m i= q K £ k sin5 + (S 0 + b w )} (4) 

where 

= E k + (to/k) Bj tan a. 
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The introduction of the variable v may be viewed as a mathematical change of 
variables, and does not presuppose the existence of a physical coordinate sys- 
tem moving with velocity (co /^Sa- 
lt follows easily from (4) that the energy 

E = -^-m v 2 + (q/k) cos 0 

is a constant of the motion. 


In a cylindrical coordinate system with axis parallel to B Q , (4) has the 
following components: 


d v 


II 


d t 


(q/m) {E^' cos a sin 9 - v^ [B_ cos (9 - ip) - B + cos (9 + i/>)]} 


^Jj 

dt 


(q/m)-|j^E^ sin a + V|| B_^ cos -^yE^ sina + v^ B + j cos + 


^ = (q/m) [B ly tan a cos 9 - B Q ] + ^ E^ sin a + v (| B_^ sin (9 - \p) 


+ f-j sin a + v (| B + j sin (9 + i//) 


(5c) 


where 

V 1 = ( V x + V y) 1/2 

'P = tan " 1 (v y /v x ) 
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B = — (B. ± B, ) 

± 2 lx ly 7 

and B. . = B„, * e .. In addition, we have the equation 

lj -V, W A 

— = k-v = k V sin vp + k v (6) 

j ^ ^ y ± z ii 

For nonresonant particles, these equations can be solved to first order 
in the wave amplitude by integrating the right hand side of equations (5a-c), (6) 
along the unperturbed trajectories, defined by 

v lf 0) = v llo> V 1 (0) =V 10 ’ 

'/' (0) = % - n 0 t , 

5<0) = e 0 + ( k y v l0 /n 0 )[cos ( 4 > q - n 0 t) - cos i/» 0 ] + k z v |10 t , 

where = qB Q /m and v ||Q , v, Q , are constants. The first order 

correction terms calculated in this manner oscillate rapidly in time and remain 
small. 

Orbit expressions calculated in this way for resonant particles contain 
slowly varying, growing terms which cause the perturbation expansion to 
break down after a finite time, as the actual particle trajectories drift farther 
and farther away from the unperturbed trajectories. This breakdown may be 
avoided by defining a new zero order trajectory. One replaces the constants 
v . . . & 0 by parameters which vary slowly in time in such a way that the 
new zero order trajectory remains close to the actual trajectory. 10 
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where 


V | ( | 0 > = *|| ( T >» v i°' ) =v i (r) 

'/' (0 ' ) -P o CO - 1 * 

0(0 ') = 0 (T) + [k y Vj_ (T)/n 0 ] COS bp (r) - n 0 t] 

- ( k y V 10 /n o) COS ^0 

+ k 2 v„ (r) t , 

T = fc » < < L 


For example, if we substitute 

V = v(0 (r, t) + v ( *) (r, t) + . . 

0 = 0(° ) (r, t) + 0( 15 (r, t) + — 

in (5a), discard terms of second or higher order in the wave amplitude, and 
Fourier expand the right hand side with respect to t, we obtain 


dv m 

d t 


“ 1 7, 


dv H ( T ) 

0 dr 


±00 


cos a J n (Uj.) 


+ \ [B_ J n+1 C^) + B + J N _j (u ± )] } sin [(k z v„ - n n Q ) t + yj (7) 


where v Ml = e • v ( 1 >, 

hi ~ ^ 7 


Z 'V/ 


u, (t) = k y v i (r)/n n , 


Y n ( T ) =5 (r) + n 


* V > 77 

P ( T ) + -j 


- ( k y v 10 /fi 0 ) cos P 0 ’ 


( 8 ) 


(9) 
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and J n is the Bessel function of order n. When the resonance condition 

k z V|, - N f2 0 = k, V„ - O) + N n 0 ~ 0 

is satisfied, the acceleration dv^ /dt has a slowly varying component (the n = 

N term in the series) which produces secularities in Vyj unless we choose 
v 1( (t) such that 

d v„ (r) _ 

V 0 — = — {Ek cos a J (u, ) 

u d r m K w - 

tv x (B. J„ n (ijl.B,],., )] } s i n $> N (r) ( 10 ) 

where 

( T ) = ( k z - n n 0 ) t + y N (t). (ii) 

Given (10), one can integrate (7) with respect to t, treating the slow variable r 
as a constant, to obtain 


± cu 

II {E k COS Jn + Yl [B - J n +l ^1> +B + h-1 C^)]} 


n^N 


{cos y n - cos [(k z v H - n n 0 ) t + yj} 

( k * v„ - n n n ) 


(12) 


Equations for the remaining slow variables, 


dv, 


V, 


o d 


r m 


■jv n [B_ J n+1 (u^ +B + J NM (u L )) 

\K sina [ T n+ i (u x ) +.T n -i Sin \ 


(13) 
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V, 


d 4> _ q 


0 d r m ] ^ 


[B_ J N+1 (u x ) -B + J NM (up] 


(2^)^ Sin a CJn_1 ^ " Jn+1 ^ l)] 


B ly tan a J N (up V cos 0 N , 


( 14 ) 


d e 

^Od7 = ^o tk 


dV ll q pz -II _ 

■tt + m B ^ tan Hit" j n ^i) 


+ | [ Vi (“1) + Vi ^i> ] | cos K 


(15) 


and for the first order corrections to the velocity components v A and </» , 



tB- J„+i («.) 


+ Jn-l 


(up] 


+ — E.' s i n a 
2 k 


[J 


n+ 1 


( S P + J n - j (^P] 


''cos [(k z v B - n n 0 ) t + y n ] - cos yj 

( k z Vy - n n Q ) 


(16) 


*1 




(up -B + J nM 



1 ^y^ E k sin a [J n+1 (u A ) - (up] 
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(17) 


+ B ly tan 


a Jn ( u 1 )|- 


fs in [(k z V|| - n Q 0 ) t + y n ] - sin yj 
( k z V|| - n n 0 ) 


are obtained in the same way. 


Reduction of Phase Space 

Only three different dependent variables (v^ , , and <£ N ) appear in the right 

hand sides of (10), (13), (14), and (15). With the aid of (9) and (11), equations 
(14) and (15) may be combined to yield a single equation for d0 N /dr . Then, after 
using the recursion formulas for the Bessel functions and letting r-17 t, one 
obtains the following closed set of equations: 


IT = mi F * COSaJ N (Sf l) 


+ V, 


B 


lx 


n 


k V. 
y 1 




■sin 4 > (18) 


d T 

d t 


m I Vkv. 


•Tn 


+ v„ 


/ n a\ 

y 1 


sin 0 


(19) 






ly 


k y V|| / N 2 n 0 2 


n 


0 


k 2 V , 2 

y 1 
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( 20 ) 


+ N I 


( fl. k N 

k z/ 


J N (“i)r cos ^N 


where J'(u. ) = dJ /du. , u. = k v, /Q . Once v M , v, , and $> are known, 6 
n' 1 ' Nil y 1 o II 1 N 

and may be obtained immediately from (14) and (15) by quadrature. 


Exact equations for a particle in resonance with an electrostatic wave 
propagating along B 0 may be recovered from (18)-(20) by setting (B ix , B iy , a » 
k y , N) = 0. Similarly, (18)-(20) reduce to exact equations for a particle in 
resonance with a parallel propagating, right circularly polarized transverse 
wave, in the limit (E/, a , k ) = 0, B, = B, , N = 1. In these limits there are 

' k y' lx ly 

no rapidly oscillating fields acting on a resonant particle, hence the first order 
corrections vanish. 


Solution of the Slow Variable Equations 
It follows from (18) and (19) that 


_1 ^ = (q/k) E ‘ (kz ^ _ N n 0 ) J N (U ± ) Sin 4> u 

where v 2 = v" 2 + Vj 2 . But, from (20), 

d 0 

(k v.| - N Q.) = + 0 [(amplitude) 1 ] 

z ii u d t 

and, from (19), 


(21) 


( 22 ) 


-d- [(q/k) J N (up] =0 [(amplitude) 2 ] 
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Hence 


(q/k) (k z v || - N n„) J N (u,) sin$ N 

= t- (q/k) E^ J N (u x ) cos $ N ] + 0 [(amp. ) 2 ] (23) 

If we combine (21) and (23) and define 


E N 3 ^ + (q/k ^ ^ (“l.) C ° S ^ 24) 

we find that dE N /dt ~ O [(amp.) 2 ]; thus E N is a constant, to the present order 
of approximation. E N is the energy associated with the slow or averaged motion. 
If one begins by differentiating m (v^ - N fi 0 /k 2 ) 2 /2 instead of m v 2 / 2, and 
then follows the same procedure used above, a second constant is obtained: 

" J m (''ll ' N Q 0 /k z )2 + (d/ k 2> COS a J N (“]) 


+ V, 


B 


lx 


n a 


k v, 
y 1 


-T n (V +B 1vTn(Si) 


cos 


(25) 


E^ and W n become exact constants of motion in the limit of parallel propagation, 
with N = 0 for an electrostatic wave and N = ±1 for a circularly polarized trans- 
verse wave. 


With the aid of (25) and (22) the resonance condition may be given somewhat 
more precisely than before: a particle will be termed resonant if the kinetic 
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energy m (v - Nf^/k z ) 2 /2, which is = m k~ 2 (d $ N /dt) 2 /2 by (22). is com- 
parable to or less than the amplitude of the effective potential defined by the 
second term on the right hand side of (25). Thus the resonance condition is 


where 


k z v.| - 0 (7?) 


V (v.) - 


(q k z /m fi 




cos a J N (u x ) 


(26a) 


+ V 1 


B 


lx 


n n. \ -n 

<“l> +B ivTn (\)\\ 


1/2 


(26b) 


The appropriateness of this definition can be seen by solving (25) for v^ and 
subtracting initial values: 


V|| (t) - v |10 = [2 W N /m - 2 a (r, fi Q /k z ) 2 cos £„] 172 |* (27) 

r \ J > 

where n , is the sign of the coefficient of cos 0 N in (25). If (26) is satisfied 

T) 

the quantity in square brackets is of first order, hence 

V|, (t) - v ||0 <0 [(amp.) 172 ] . * 28) 

When (26) is not satisfied W N is much greater than the term proportional to 
t) 2 in (27), and we can expand the square root, to obtain 

v,| (t) - v ||0 <0 [(amp.) 1 ]. 

Thus (26) correctly gives the condition for an enhanced or resonant interaction. 
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It follows immediately from (24) that 


v 2 (t) - v Q 2 <° [(amp.) 1 ], (29) 

and it follows from (28) and (29) that 


^ (t) - v lQ < 0 [(amp. ) 1/2 ] . (30) 

To lowest order in the square root of the wave amplitude one can ignore the 
small variation with time of ^ (t) in terms which are already of first order, in 


(24), (25), and (18). Thus 


m v ■ 


(q/k) E k ' J N (u i0 ) COS 


(31) 


W n = o m - N 0)/k z ) 2 + m a (rj fi 0 / k ) 2 cos 0 


(32) 


d t 



where % = V (v l0 ). 


(33) 


The path of a particle in the two dimensional phase space with coordinates 
$> N , V|| can be determined by solving (32) for 0 N and plotting versus v . 

A typical set of such paths, for a fixed value of v, 0 and several different values 
of W N , is illustrated in Fig. 1. Note that the angle is bounded for particles 
following the closed orbits. Such particles are said to be trapped. 
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Eq. (32) may also be used to eliminate $> N from (33). The resulting equation 
may be written in the form 


^=±[(l-Y 2 )(Y J + r*-l)]>/’ < 34 > 

d s 

where 

v = ( r /2 V (k, v H /n 0 - N) 
d s = - (T7 0 n 0 /r) d t 

r 2 = 2 (1 + k 2 W N /m 17 2 f) 2 )" 1 ( 35 ) 

Eq. (34) is the equation satisfied by the Jacobian elliptic function =f dn (s, r). Thus 


v„ (t) = Nn 0 /k z +(2 t] 0 n 0 /r k z ) dn [(^ fi 0 /r) (t + - t), r] (36) 

The function dn (s, r) is periodic in s. The period of the motion, for untrapped 
particles, is 

T„ = 2 r K (r)/, 0 (^) Q, < 3 ’> 

where K(r) is the complete elliptic integral of the first kind. For untrapped 
particles with velocities greater (less) than Nf^/k z one chooses the lower 
(upper) sign in (36). The phase constant t + is any time at which the quantity 
|v|| - Nf^/kz | has its maximum value. The parameter r classifies the particle 

orbits according to type, r is greater than one for trapped particles, with r - 00 
for particles in the center of the nested set of closed orbits in Fig. 1. r has 
the value unity for particles on the separatix and is less than one for untrapped 
particles, with r < O [(amp) 172 ] for nonresonant particles. 13 
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Strictly speaking, (36) is valid for both trapped and untrapped particles, since 
dn(s, r) is defined by analytic continuation when r > 1. It is convenient, however, 
to reserve the use of (36) in its present form for untrapped particles (r < 1), and 
to use the identity dn(s, r) = cn(rs, l/r) when r 2 1. Thus 

''ii (*) = N n„/k, + (2 77„ n o /r k 2 ) cn (t, - t), l/r] < 38 ) 

t k = 4K(i/o/^ (v 10 >n 0 (39 > 

for trapped particles. We have omitted the upper sign on the last term in (36), 
in deriving (38). Adding one half period to t + changes the sign of the elliptic 
function cn, hence no generality is lost through this omission. 


After solving (32) for $> N and using (35) and (36), one has 

(t) = cos -1 f(2 cr 2 /r 2 ) {1 - dn 2 [(17 Q. /r) (t - t + ), r]} - a \ ( 40 ) 

V V T) / 

The time dependence of is found by solving (31) for v^ 2 , subtracting initial 
values, and using (36) and (40). After taking and expanding the square root of the 
resulting equation, one has 


Vj. (0 - v iQ ± (2 N v 0 H 2 /r k 2 v i0 ) dn fi 0 /r) (t - t + ), r] 

to lowest order in (amp.) 1/2 . Again, one can use the identity dn(s, r) = 
for trapped particles. 


cn (rs, l/r) 
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In summary, the major effect of the resonance on the motion of single 
particles is to induce the particles to execute large amplitude amp. 1/2 ), long 
period (~ amp. 1/2 ) oscillations about their zero order trajectories. Super- 
imposed on these large motions are rapid, small amplitude oscillations. 

III. SATURATION OF RESONANCE DAMPING 
The physical mechanism underlying the Landau damping of electrostatic 

14 

waves was described by Dawson in 1961. Dawson demonstrated that the 
damping was due to the presence of resonant particles which move with velocities 
very close to the wave velocity. These particles exchange energy with the 
wave on a time scale long compared with the period 2tt/co. If the resonant 
particles are initially distributed or bunched in phase space in such a way that 
more particles gain than lose energy, then the wave damps, by conservation 
of energy. The kinetic energy of a resonant particle in a sinusoidal potential 
oscillates in time, but the period of oscillation approaches infinity if the wave 
amplitude is infinitesimal. Thus a very small amplitude wave damps monotonically. 

If the wave amplitude is large the periodic nature of the particle motion 
can not be neglected. This situation has been treated by O’Neil. 1 After a time 
of the order of a half period of a typical resonant particle's motion, and if the 
wave amplitude has not changed too much in the meantime, most of the particles 
which were initially gaining energy begin to lose energy, and vice versa. The 
damping coefficient then changes sign, and the wave amplitude begins to oscillate. 
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The bunching of resonant particles in phase space, required for a net exchange 
of energy between wave and plasma, cannot be maintained indefinitely. The fre- 
quency of oscillation of a particle in a potential well of the wave depends on 
the particle's energy, and the spectrum of oscillation frequencies is broad. Thus, 
phase mixing of resonant particles eventually causes the damping coefficient to 
vanish entirely. In the absence of instabilities one is left with a wave of constant 
amplitude. 

The cyclotron damping of a finite amplitude whistler wave propagating parallel 
to a background magnetic field has been shown to saturate in the same way. 2 

In the preceding section we have shown that the motion of particles in 
resonance with an obliquely propagating wave exhibits, to lowest order in (amp.) , 
all of the features relevant to the nonlinear damping process described above. 

That is: 

1. resonant particles exchange energy with the wave on a time scale long 
compared with 2tt/u>\ 

2. the kinetic energy of a resonant particle in an oblique wave is periodic 
in time; 

3. the spectrum of resonant particle oscillation frequencies is broad. 

One therefore expects that the damping of a finite amplitude oblique wave should 
also vanish after some time, to lowest order. 
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In this section we shall describe the resonance damping of an obliquely 
propagating wave in the nonlinear regime. Our purpose is not to provide a 
fully detailed picture of the damping process, but rather to determine the time 
dependence of the nonlinear damping coefficient in a simple, somewhat heuristic 
way. Specifically, we shall assume that the damping coefficient obtainable from 
the linearized theory, 



(41) 


N 

is known, and describe the nonlinear coefficient y(t) in terms of functions h N (t), 
defined by the relation 


| y(t) dt = 22 


(42) 


The index N numbers the contributions from the various resonances. 


Following the method used in Ref's. 1 and 2, we begin with the energy conser- 
vation equation. If the amplitude parameters E k , B. have a time dependence of 
the form 


A (t) = Aq exp 



7 (t) d t , 


then the energy equation in the lab frame is 


2Qi 7 (t) +P nr (t) + 


F [ P r(t)] N = 0 


(43) 
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where 


Qi 


2tt e Q 



(B j + B l) 


The quantity 2Q y is the time rate of change of the wave energy per wavelength. 
P nr and 2 (P r ) N represent the rate of energy absorption per wavelength by non- 
resonant and resonant particles, respectively. 

The time scale on which the wave amplitude varies is assumed to be much 
longer than the period of oscillation of a nonresonant particle, hence adiabatic 
invariants exist which relate the constants of motion of a nonresonant particle 
at any given time to the instantaneous value of the wave amplitude. If follows 
from this that the total kinetic energy per wavelength of nonresonant particles 
depends on time only through the wave amplitude. Hence 


P nr (*■ E ->„,=Q 2 ?(t) < 44 > 

where Q 2 is a function of the wave amplitude and various constants. 

After using (44) in (43), the latter equation may be solved for y (t). 

7 (t) = - L +Q,)1 <i p , <t)] N /(P,U (45 > 

N 

where (P r ) LN is the value of [ P r (t)J N obtainable from the linearized theory. If 
we recall that linearized theory provides an adequate description of the behavior 
of nonresonant particles, in the present approximation, it follows from a com- 
parison of (41), (42), and (45) that 
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(46) 


+ ^2^ ^LN 

< t >V< P A.N = ^ h »< t >' 

The rate of energy absorption by resonant particles is given by the relation 


[p < (t >k = f- t 


' 5 " 


IK d3 v I d e 

V.. J o 


2 u> v„ 


(v 2 + 


a>- 


k 2 


f (#> ▼. t) 


(47) 


The quantity in parenthesis in the integrand above is the square of the lab frame 
velocity, f is the distribution function, and V N is the volume in velocity space 
associated with the Nth resonance. To estimate the relative magnitude of the 
various terms in (47), it is helpful to write [ P r (t)] N in the alternative form 


fe, (*)]* = 


1 . 

2 



[v 2 (ff 0 . 


Z O' 


t) 


+ 2wv H (6 0 > Z o’ t) /k z +" 2 /k 2 ] f (0 Q , v 0> 0) 

The constant m 2 /k 2 vanishes under the time derivative and makes no contribu- 

Z 

tion. If saturation occurs before the wave amplitude has changed very much the 
approximate trajectories developed in Sec. II may be used to express v 2 (<9 Q , v , 
t) and V|| (9 q , v , t). v 2 and Vj| may be separated into fast and slow parts. As 
far as damping is concerned, the important effect of the resonance is to produce 
the large, slow deviations from the zero order particle trajectories described 
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by the slow variables. Thus we consider only the slow parts. It follows from 
(18) that dv (| /dt is a first order quantity, while (21) and (36) indicate that dv 2 /dt 
is of the order (amp.) 3/2 . (After integration the term is of order (amp.) 2 , 
while the v 2 term is of order (amp.) 5/2 .) Thus, to lowest order, 



We shall assume that f is Maxwellian when the wave is "turned on" at t = 0. 
The small deviation from thermal equilibrium required to support the wave 
initially will be neglected. In general, this initial perturbation phase mixes away 
rapidly and makes no significant contribution to the damping. In a collisionless 
plasma the phase space density in the vicinity of a given particle is conserved, 
hence 


f (^> v, t) = A exp {- [v 2 (0, v, t) + 2wv ||0 (6\ v, t)/k z 


+ co 2 /k 2 l/V 2 } 


where v 2 , V|| 0 are the values of v 2 , Vy at t = 0 for a particle at 0, v at time t, 
and A is a normalization constant. The right hand side of (49) will be expanded 
in the wave amplitude. It follows from (28) and (29) that the lowest order time 
dependent term in the expansion will be of order (amp.) 1/2 and will come from 
v ||0 ( 6 , y, t). To this order, v 2 ~ v 2 . Thus, after setting v |(0 equal to NQ 0 /k z 

plus the small quantity (v - Nfi o /k z ) and expanding, we have 
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f (6, v, t) = A exp {- [v 2 + (a? + 2NwQ 0 )/k 2 ] /V 2 } 

x [1 - (2^/k z V 2 ) (v ||0 - Nfi 0 /k z ) -• • •] (50) 

The integral in (48) can be evaluated if one additional simplifying assumption 
is made. We shall assume that the perpendicular velocity of a typical resonant 
particle is small compared with /k y , and neglect terms of order (k y vj_ /fi Q J 2 = 
u 2 . Thus, for N = 0, (26b) gives 

7 ?o s ^( v l 0 ) = ta k z E k cos a/mQ 2 | 1/2 + 0 (u 2 q ). ( 51 ) 

WhenN = ±1, tj 0 is proportional to v^ 1 / 2 . For a resonant particle = [v 2 - 
(NQp/kz) 2 ] 1/2 + O [(amp.) 1/2 ] . Thus, when N = ±1 

r 1 1/2 

V 0 = (q/mfi 0 ) j (B ly ± B lx ) + k y E' cos a/2fi Q 

x |(k z v/Q 0 ) 2 - l| 1/4 + 0 [(amp. ) 3/2 , u£ 0 ], (52) 

When | N | > 2, v 0 - 0 (u^ 0 ). It follows from this last relationship and (36), 

(45), (48), and (50) that 

\ ' i dh N 

7<t) = 2_, y u> -57 + 0(<u l» > ) < 53) 

N=~l 

When N = ±1, the oscillatory part of the motion described by (35), (36), (38), 
and (52) is identical to the approximate parallel oscillatory motion of a particle 
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in resonance with a whistler wave propagating along the magnetic field, if the 
wave number and amplitude of the whistler are given by 

k = k 

eff. z 

B. = — (B. + B, ) + k E' cos a/2fl A . 

[See Eq.'s. (21a), (21b), and (22) of Ref. 2.] In addition, (48) and (50) are 
identical in form to the corresponding equations [(25), (26') and the first eq. 
in Appendix B of Ref. 2.] for the parallel propagating whistler. The rate of 
energy absorption by particles in resonance with a parallel propagating whistler 
is calculated explicitly in Ref. 2. We can therefore obtain expressions for 
[P. (t)] by making appropriate generalizations of the constants in the results 
of this earlier work. 

Similarly, it is easily shown that (35), (36), (38), and (51) make the evaluation 
of (48) in the case N = 0 mathematically equivalent to the calculation of P (t) for 
a finite amplitude electrostatic wave. We can obtain [P (t)] Q by making approp- 

l 

riate generalizations of the constants in O'Neil's result. 

Finally, it is demonstrated in Ref's. 1 and 2 that, for N = 0, ±1, 

( P r) L N =limit (54) 

t/T N*° 

where r is the period of oscillation (slow) of a typical resonant particle. 
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After following the procedure for finding [P. (t)] Q ±1 outlined above and 


using (46) and (54), we have 


r 


'« (t) = t » tL t J dr 

-a 


277 1 - COS 


< 


77n t 


L r K ( r ) T 0 j 


r 4 K(r) (1 + q 2n ) (1 + q~ 2n ) 


2n > 


2 tft l 1 - COS 


(2n + 1) 77 1 


L 2K(r) r Q J, 


K(r) (1 + q 2n+1 ) (1 + q" 2 "" 1 ) 


where q = exp {- ttK ( [1 - r 2 ] 1/2 )/K(r)} and 


(55) 


= | (qk/m) E' cos 2 a|~ 1/2 ; 


and 


h ±1 (t) = 2 9 tt 1/2 - 1^1 T±1 


1 


-X 


x| dr 

\ n=l 


£ 


1 - exp {-[(n77-/2C x C 0 rK)(t/T ±1 )] 2 } cos [(nw/rK) (t/r ±1 )] 
r 4 K(r) (1 + q 2n ) (1 + q" 2n ) 


00 

L 

n- 0 


[1 - exp (- { [(2n + 1) tt / 4Cj £ q K] (t/r ±1 )} 2 ) cos { [(2n + 1) 7 t/ 2K] (t/r ±1 )}] r 


K(r) (1 + q 2n+1 ) (1 + q“ 2n_1 ) 


(56) 
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where 


±1 


( qk z 


a/™> [i 


1 - 1/2 


( B i y ± B lx ) + k E' cos a/2Q ( 


and £ 0 ~1.1, C 0 ~ .41, C t % 3.1. 

When these results are substituted into (53) our estimate of the damping 
coefficient is complete. The time variation of h Q (t) has been described in the 
introduction to this section. The behavior of the functions h x (t) and h_ x (t) is 
similar, and is sketched in Fig.'s 3 and 4 of Ref. 2. 


The asymptotic values h Q ±l (oo) have been calculated in the papers quoted 
above. One finds that 


f 


y(t) dt ~ 2 


+1 

L 

N=-l 


^LN T N 


In deriving (55) and (56) we have assumed that saturation occurs before the 

J '.ao 

ydt « 1 

0 

In a multi -component plasma, these results represent the contribution of 
any one species to the total damping coefficient. 


IV. SUMMARY 

We have found approximate solutions, of relatively simple form, to the equa 
tions of motion for a particle in resonance with a small finite amplitude wave 
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propagating at an angle to a magnetic field in a plasma. These solutions can be 
applied to the study of several important effects which involve resonant particles. 

One such effect, resonance damping, has been considered here. We have 
shown that the resonance damping coefficient for a finite oblique wave vanishes 
after a few bounce times, leaving, in the absence of instabilities, an undamped 
wave. 
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